CRISPR/Cas9-mediated disruption of CjACOS5 confers no-pollen formation on sugi trees (Cryptomeria japonica D. Don)

Sugi (Cryptomeria japonica D. Don) is an economically important coniferous tree in Japan. However, abundant sugi pollen grains are dispersed and transported by the wind each spring and cause a severe pollen allergy syndrome (Japanese cedar pollinosis). The use of pollen-free sugi that cannot produce pollen has been thought as a countermeasure to Japanese cedar pollinosis. The sugi CjACOS5 gene is an ortholog of Arabidopsis ACOS5 and rice OsACOS12, which encode an acyl-CoA synthetase that is involved in the synthesis of sporopollenin in pollen walls. To generate pollen-free sugi, we mutated CjACOS5 using the CRISPR/Cas9 system. As a result of sugi transformation mediated by Agrobacterium tumefaciens harboring the CjACOS5-targeted CRISPR/Cas9 vector, 1 bp-deleted homo biallelic mutant lines were obtained. Chimeric mutant lines harboring both mutant and wild-type CjACOS5 genes were also generated. The homo biallelic mutant lines had no-pollen in male strobili, whereas chimeric mutant lines had male strobili with or without pollen grains. Our results suggest that CjACOS5 is essential for the production of pollen in sugi and that its disruption is useful for the generation of pollen-free sugi. In addition to conventional transgenic technology, genome editing technology, including CRISPR/Cas9, can confer new traits on sugi.

www.nature.com/scientificreports/ induce type II immunity. Japanese cedar pollinosis was first reported in 1964 in Nikko, Japan 5 . Thereafter, several nationwide surveys showed that the number of patients with Japanese cedar pollinosis had been increasing: the estimated prevalence was 11.7% in 1998 6 , 13.1% in 2001 7 , 26.5% in 2008 8 , and 38.8% in 2019 9 . The reduction of the amount of sugi pollen has been thought as a countermeasure to Japanese cedar pollinosis. For this purpose, wild varieties with fewer male strobili and wild male-sterile mutants without pollen grains have been used for breeding. Male-sterile mutants were discovered in many plant species and the genetic patterns of male sterility are divided into genic male sterility (GMS) and cytoplasmic male sterility 10,11 . GMS is also called nuclear male sterility or Mendelian sterility because it depends exclusively on the nuclear genome and exhibits Mendelian inheritance. The development of genetic analyses has served to identify several key genes involved in GMS. For example, a transposon insertion mutant of the Arabidopsis thaliana acyl-CoA synthetase 5 (ACOS5) gene does not produce mature pollen in anthers 12 . The recombinant ACOS5 protein catalyzes oleic acid CoA-ester formation in vitro, and thus probably participates in the biosynthesis of sporopollenin, which is a constituent of the exine of pollen grains. Rice (Oryza sativa) carries the OsACOS12 gene, which is an ortholog of ACOS5, and the osacos12 mutant also shows a phenotype of absence of mature pollen 13,14 .
A male-sterile sugi without pollen was first identified in Japan in 1992 15,16 ; subsequently, various male-sterile mutant trees were discovered and analyzed [17][18][19] . To identify the genes that are responsible for male sterility in sugi, genetic marker analyses have been performed, leading to the identification of four recessive male sterilitylinked independent loci (ms1, ms2, ms3 and ms4) [20][21][22][23][24] . Recently, one gene (MS1) at the ms1 locus was identified as a causative gene for male sterility 25 . MS1 encodes a lipid transfer protein and is expressed in male strobili, but its biochemical properties or physiological roles remain unknown. The appearance rate of pollen-free sugi trees is approximately 0.02% in a seed orchard in Japan 26 , and these trees are extremely rare in the field. In addition, the discovered mutant trees should be crossed to sugi plus trees with excellent properties for the generation of pollen-free plus trees, because most of the original male-sterile mutant trees exhibit slow growth or curved trunks. Therefore, the development of male-sterilization technology is valuable and needed for the generation of pollen-free sugi plus trees.
In this study, we aimed to generate pollen-free sugi trees using the CRISPR/Cas9 system. Genome editing technology including the CRISPR/Cas9 system has been widely utilized in plants and animals for the elucidation of gene function and for the breeding of new varieties. Because of the simplicity of its experimental procedures, the CRISPR/Cas9 system has been applied not only to herbaceous plants and crops, but also to woody plants such as sweet orange trees 27 , apple trees 28 , grape vines 29 , poplar trees 30,31 , sugi trees 32 , radiata pine 33 , and white spruce 34 . Although the use of genetically modified organisms (GMOs) including transgenic trees have been regulated by law in Japan, the products by the genome editing technology are less restricted by law, if they have no DNA from other species. Recently, γ-aminobutyric acid (GABA)-enriched tomato which had been generated by CRISPR/Cas9 started to be commercially sold without the regulation against GMOs 35 . The genome editing technology has the potential to be put the modified and improved forest trees to practical use in future. Here we report the genome editing of sugi by the CRISPR/Cas9 system, to mutate the CjACOS5 gene, which is an ortholog of Arabidopsis ACOS5 and rice OsACOS12. The regenerated sugi trees with mutated CjACOS5 genes flowered but did not produce pollen grains in male strobili. We demonstrated that the CRISPR/Cas9 system enables the generation of pollen-free sugi trees artificially, even though we had not been able to achieve this in the past.
Furthermore, we cloned cDNAs homologous to CMFL003_A04 from male strobili of sugi #13-8-2 line and sequenced. As a result, we identified two cDNAs, CjACOS5a (accession number: LC726337) and CjACOS5b (LC726338). The predicted protein encoded by CjACOS5a and CjACOS5b comprises 557 amino acid residues and their 551 residues (98.9%) are equal to each other (Fig. 1). The CjACOS5a protein has 61.9% identity to ACOS5 and 59.8% identity to OsACOS12. The CjACOS5 proteins are composed of two domains; one is an AMPdependent synthetase/ligase domain (InterPro ID: IPR000873), from Glu41 to Tyr454, and the other domain is a C-terminal domain in AMP-binding enzymes (IPR025110), from Glu463 to Lys538 (Fig. 1). Both domains also exist in ACOS5 and OsACOS12. In the AMP-dependent synthetase/ligase domain, the CjACOS5 proteins have a conserved AMP-binding site (IPR020845), from Leu197 to Lys208 (LPYSSGTTGASK) (Fig. 1). This conserved AMP-binding site also appears in long-chain fatty acid-CoA ligases and 4-coumarate-CoA ligases, according to the InterPro web site. OsACOS12 also carries the conserved AMP-binding site (197-LPYSSGTTGVSK-208), whereas ACOS5 carries a similar amino acid sequence (186-LPFSSGTTGLQK-197), which was not determined to be IPR020845. As the result of these genetic analyses, CjACOS5a and CjACOS5b were thought to be sugi orthologs of ACOS5 and OsACOS12 and encode an acyl-CoA synthetase. The expression of the CjACOS5a gene in male strobili was about 280-fold higher than that detected in stems, whereas that in the leaves and female strobili was not significantly different from that in stems ( Supplementary Fig. S1). Similarly, CjACOS5b was expressed most highly in male strobili, and their mRNA level was about 150-fold higher than that in stems. Accordingly,

CRISPR/Cas9 vector construction and sugi transformation. To induce loss of function of CjACOS5
in sugi, we planned the target sites for DNA breaks by the CRISPR/Cas9 system. The recruitment of the complex of guide RNA and Streptococcus pyogenes Cas9 requires 5′-NGG-3′ as a protospacer adjacent motif (PAM) and the upstream target sequence. The AGG sequence located downstream of the start codon of CjACOS5 was choose www.nature.com/scientificreports/ as a PAM (Fig. 2a). If the CjACOS5 gene is broken at the predicted site, at 3 base pairs upstream of the PAM, short CjACOS5 proteins with about 34 or 24 amino acid residues would probably be synthesized by a flame shift. Two allelic CjACOS5 genes, CjACOS5a and CjACOS5b, were distinguished by a single nucleotide polymorphism (T or G) located upstream of the start codon (Fig. 2a). A CRISPR/Cas9 vector was constructed for transformation of sugi. The original pZK_gYSA_FFCas9 vector was generated for dicots and is composed of the A. thaliana U6 promoter-driven guide RNA, the parsley ubiquitin promoter-driven Cas9, and the cauliflower mosaic virus (CaMV) 35S promoter-driven NPTII gene 37 . We modified pZK_gYSA_FFCas9 to generate the pBFGE1 vector, in which the Zea mays ubiquitin promoter was inserted between the CaMV 35S promoter and the NPTII gene, to further express the NPTII gene (Fig. 2b). The ubiquitin promoter has been reported to work well in sugi 38 . The 20-bp target sequence of CjACOS5 was integrated into the guide RNA gene of the pBFGE1 vector, to break the CjACOS5 genes (Fig. 2b).
The sugi embryogenic cell line #13-8-2 was infected with Agrobacterium tumefaciens GV3101 harboring the constructed pBFGE1 vector to generate transgenic sugi trees. After infection, the cells were cultivated in selection medium containing kanamycin and formed calli harboring the vector. We picked six proliferated calli and investigated whether they had the NPTII gene by PCR. Among the six calli, three kanamycin-resistant calli lines were obtained and designated as GE#1, GE#2 and GE#3 ( Supplementary Fig. S2a). These kanamycinresistant calli and the non-transgenic #13-8-2 calli were transferred to kanamycin-free maturation medium, and then all of the calli formed somatic embryos in about 2 months. Subsequently, the somatic embryos were moved onto germination media, germinated, and grown. The regenerated transgenic plantlets showed normal shape with needles and roots as same as regenerated non-transgenic plantlets. The rooted transgenic sugi plantlets were planted in a pot and grown in a closed glass greenhouse (phytotron). All sugi plantlets regenerated from Figure 2. Design of the binary vector for the CRISPR/Cas9 system. (a) CjACOS5 cDNA sequence neighboring the DNA break site for the CRISPR/Cas9 system. The target sequence and PAM are shown in red and underlined blue, respectively. The arrow head is located at a predicted DNA break site. The reversetyped K is a single nucleotide polymorphism and indicates thymine in the CjACOS5a gene or guanine in the CjACOS5b gene. Amino acid residues are indicated under the cDNA sequence. (b) Schematic representation of the constructed CRISPR/Cas9 binary vector. This vector was named pBFGE1 and was derived from pZK_ gYSA_FFCas9 37  www.nature.com/scientificreports/ the GE#1, GE#2 and GE#3 calli possessed the NPTII gene in their leaves, whereas non-transgenic sugi plantlets obtained from the #13-8-2 cells did not ( Supplementary Fig. S2b).
Induced mutation of CjACOS5. To confirm the mutation of CjACOS5 in the transgenic trees regenerated from the GE#1, GE#2 and GE#3 calli, genomic DNA was isolated from the leaves of each transgenic sugi plantlet, then DNA fragments including the target site of CjACOS5 were amplified using PCR and cloned into a plasmid vector in E. coli. The plasmids were individually sequenced and analyzed. As a result, several deletion or insertion mutants were detected in the target DNA sequence of CjACOS5a and CjACOS5b in the transgenic sugi trees, whereas non-transgenic sugi trees that were regenerated from embryogenic cells had no mutation in the two genes (Supplementary Table S2). The two GE#1 callus-derived transgenic sugi lines, GE#1-17 and GE#1-18, were chimeric plants (Fig. 3). The GE#1-17 sugi tree carried a 3-bp deletion and a 1-bp deletion of CjACOS5a and two types of 1-bp deletion of CjACOS5b. GE#1-18 contained wild-type and 1 bp-deleted CjACOS5a and a 1-bp deletion of CjACOS5b. In contrast, all eight transgenic sugi lines derived from GE#2 calli (GE#2-1 to GE#2-22) had the same biallelic mutation with a 1-bp deletion at the same position of CjACOS5a and CjACOS5b (Fig. 3 and Supplementary  Table S2). These deletions of a single base pair (G:C) suggested the production of a truncated and mutated CjACOS5 protein of 34 amino acid residues (MADQNAVREKRRSFSAAPFQLFWFLRTLVYRNSF). In this mutant CjACOS5 protein, the nine underlined-N-terminal amino acid residues alone are in accordance with the original CjACOS5 protein. Conversely, GE#3-derived transgenic trees include several types of mutation. Five transgenic sugi lines (GE#3-7, GE#3-8, GE#3-9, GE#3-11, and GE#3-12) were chimeric plants including wild-type and deleted CjACOS5 genes ( Fig. 3 and Supplementary Table S2). There were not only 1-bp deletions in those plants, but also a 45-bp deletion in GE#3-9, and a new DNA sequence comprising a 17-bp deletion and a 24-bp insertion in GE#3-12. GE#3-10 alone did not contain the wild-type CjACOS5a and CjACOS5b genes; thus, it was a biallelic mutant, similar to the GE#2 lines (Fig. 3).
We investigated off-target mutagenesis of sugi genes which had similar sequences to the CjACOS5 target sequence and the PAM. Because the sugi genome is huge and has never been clearly elucidated yet, three genes, CJt056454, CJt031386, and CJt113843, were chosen from the sugi transcriptome assembly (CJ3006NRE) 39 using the blastn analysis. In conclusion, the PCR products of the three genes from the investigated transgenic sugi lines, GE#1-17, GE#1-18, GE#2-1, and GE#3-7 had only normal sequences and showed no off-target mutagenesis ( Supplementary Fig. S7). These results confirmed that the CRISPR/Cas9 system worked in sugi cells and induced mutations of the CjACOS5a and CjACOS5b genes.  Fig. S3a). The shape of the transgenic sugi trees from GE#2 lines and GE#3 lines appeared to be the same as that of non-transgenic sugi trees ( Supplementary Fig. S3b, e-i). The initial height of trees from GE#2 lines and GE#3 lines was similar to that of non-transgenic sugi trees at 6 months, whereas GE#1-17 and GE#1-18 trees were significantly shorter than non-transgenic trees ( Supplementary Fig. S3c, d, j). Unfortunately, GE#1-17 and GE#1-18 trees withered and died in the winter of 2017, probably because of severe growth suppression.
Flowering of sugi can be induced by spraying a gibberellin A3 (GA 3 ) solution onto the twigs with leaves 40 . At the end of July/beginning of August of 2017, the young transgenic sugi trees with the mutation of CjACOS5 and non-transgenic sugi trees were treated twice with GA 3 and cultivated in the phytotron. Visible flower buds appeared in October and developed to male strobili. In January of 2018, yellow-powdered pollen was packed in the male strobili of non-transgenic sugi trees (Fig. 4a). Spherical pollen grains in the male strobili were detected under a microscope.
In contrast, GE#2 sugi trees from GE#2-1 to GE#2-6, which were biallelic mutants of CjACOS5a and CjACOS5b, produced male strobili without pollen grains (Fig. 4b, c, and Supplementary Table S2). The male strobili were finely sliced using a razor blade in water, and the suspension was stained and observed under a microscope, however, no-pollen grains were detected and small cells were observed. Conversely, GE#3 sugi trees resulted in three types, as follows. The chimeric trees with wild-type and deleted genes, GE#3-7, GE#3-9, GE#3-11, and GE#3-12, had pollen-including male strobili and no-pollen male strobili ( Fig. 4d and Supplementary  Table S2). GE#3-8 was also a chimeric tree but had male strobili with pollen (Fig. 4e). GE#3-10, the sole biallelic mutant among the GE#3 lines, had no-pollen at all, similar to the GE#2 lines (Fig. 4f). Therefore, disruption www.nature.com/scientificreports/ of both CjACOS5a and CjACOS5b led to the loss of the ability of pollen production in male strobili. Moreover, the pollen productivity was shown to be maintained exclusively by one allele gene, CjACOS5a or CjACOS5b. A sterile male strobilus of GE#2-1 and GE#2-2 tree carried only deleted CjACOS5 genes as same as their leaves ( Supplementary Fig. S4). Conversely, a sterile male strobilus of GE#3-7 carried not only deleted CjACOS5 genes but also wild-type genes similar to a fertile male strobilus of GE#3-7. That is probably because a male strobilus had been constructed with not only mutated reproductive cells but also somatic cells without mutation. Additionally, the male strobilus of GE#3-7 possessed different types of mutations, a 2-bp deletion and a 4-bp deletion, which were not been observed in the leaves (Supplementary Fig. S4 and Fig. 3). The plausible reason is that GE#3-7 is a chimeric tree carrying the active Cas9 gene which has been inducing the mutation of the CjACOS5 genes.
To determine the reproducibility of the no-pollen flowers, the mutated sugi trees were again treated with GA 3 in the summer of 2018. They were grown in the phytotron and in a special netted-house. Biallelic mutant sugi trees, such as GE#2-1 and GE#3-10, grown in the phytotron had male strobili but did not produce pollen in 2019 ( Supplementary Fig. S5b, c, f). The chimeric mutants (GE#3-7, GE#3-9, GE#3-11, and GE#3-12) produced two types of strobili, i.e., with or without pollen (Supplementary Fig. S5d). GE#3-8 produced only fertile male strobili similar to the result of the preceding year ( Supplementary Fig. S5e). The special netted-house is bigger than the phytotron and has a glass roof and a metal mesh for ventilation. That did not have an air conditioner, so that its temperature depended on the outside temperature and sunlight but often reached around 40 °C in the summer. The sugi trees suffered no damage under high temperature conditions. In the special netted-house, the biallelic mutant trees (GE#2-21 and GE#2-22) did not produce pollen in their male strobili, whereas the non-transgenic trees had pollen ( Supplementary Fig. S5g-i). During the third experimental period, from 2019 to 2020, CjACOS5 biallelic mutant GE#2 lines also had male strobili but did not produce pollen in the special netted-house ( Supplementary Fig. S6). Accordingly, it was suggested that CRISPR/Cas9-induced mutation of CjACOS5 genes stably maintained the no-pollen phenotype over a few years in sugi trees.

Discussion
Pollen-free sugi trees have been expected to mitigate Japanese cedar pollinosis. For that reason, male-sterile pollen-free sugi trees have been the goal of time-consuming and laborious research, and are rarely found in the field or in the test forests. In the present study, we accomplished the artificial generation of pollen-free sugi using genome editing technology. Sugi CjACOS5 genes were targeted and mutated by a CRISPR/Cas9 vector via Agrobacterium-mediated transformation. The regenerated sugi trees with biallelic CjACOS5 deletions could not produce pollen in their male strobili (Fig. 4). Although small dots were observed in the suspension of no-pollen male strobili (Fig. 4), we were not sure that they were pieces of cell debris or developing microspores which were progenitor cells of pollen grains. Even if the small dots were microspores, they would be expected to be finally degraded and not to be dispersed from male strobili, because there are many reports that male sterile mutants do not release microspores outside. In the Arabidopsis acos5 mutant, microspores are lysed and completely degraded on the way of development 12 . Similarly, the rice osacos12 mutant degenerates microspores and does not have mature pollen at later stage of anther development 13 . Sugi male sterile 1 (ms1) mutants neither produce nor release pollen, because microspores are not separated in the tetrad stage and then are degenerated 16,18 . In other sugi male sterile mutants, ms2 and ms3, micropores adhere each other and collapse in anthers, and thus the microspores are not dispersed outside 17 . Accordingly, biallelic CjACOS5-mutated sugi trees are thought not to release small developing microspores outside male strobili.
Male sterility with no-pollen has been regarded as being important and useful for human health, the economy, biodiversity, and agriculture. For example, if plant species with a large amount of allergenic pollen, such as sugi, cypress, and birch trees and pasture grasses, cannot produce pollen, the number of pollen-allergic patients will decrease. As the result, the pollinosis-related health expenditures and economic loss will decrease 41 . Pollen-free plants can also avoid unnecessary gene flow, which causes problems of biodiversity in the case of invasive and GMO plants 42 . Male sterility in a variety of crops (maize, rice, wheat, etc.) has helped produce hybrid seeds more efficiently than before [43][44][45] . Using genome editing technology, it is possible to turn a plant species with pollen into pollen-free varieties, even if wild pollen-free mutants had never been found in the same plant species. Thus, the generation of male-sterile plants using the CRISPR/Cas9 system is innovative and greatly available.
The no-pollen phenotype in male strobili of the genome-edited sugi trees was thought to be caused by the mutation of CjACOS5a and CjACOS5b, as demonstrated by the results obtained for different mutation types. The GE#2 lines and the GE#3-10 line had biallelic mutation of both deleted CjACOS5a and CjACOS5b, and they showed a no-pollen phenotype in male strobili (Figs. 3 and 4). These results demonstrated that biallelic mutation of CjACOS5a and CjACOS5b was closely related to the no-pollen phenotype. GE#3-7 produced male strobili partly with pollen or without pollen. Because GE#3-7 was a chimeric tree carrying the deleted CjACOS5a gene and the normal and delated CjACOS5b gene, it was suggested that normal CjACOS5b functioned in pollen production. Male strobili without pollen and with pollen were thought to be derived from biallelic mutant cells and monoallelic mutant cells, respectively. Conversely, it was thought that CjACOS5a alone could lead to the production of pollen, because GE#3-8 had normal CjACOS5a, deleted CjACOS5b, and normal CjACOS5b and produced only male strobili with pollen. Therefore, it was concluded that the no-pollen phenotype observed in this study was caused by loss of function of both the CjACOS5a and CjACOS5b genes, and that one allelic gene, alone, either CjACOS5a or CjACOS5b, was essential to generate pollen in male strobili. If possible, to confirm that the no-pollen phenotype was not caused by off-target mutagenesis, we should elucidate the existence of off-target mutagenesis in the genome-edited sugi trees in the same way as that used for several model plants 46 . However, the off-target mutagenesis in sugi is difficult to investigate because its genomic DNA is very large (about 11 Gb) 47 and has never been adequately read. Here, we investigated the off-target mutagenesis of only three genes similar to the CjACOS5 target sequence and verified they had no mutation (Supplementary Fig. S7). www.nature.com/scientificreports/ For the development of genome editing technology in sugi, more extensive survey of off-target mutagenesis is surely needed in the future. We do not know the reason why GE#1-17 and GE#1-18 showed severe growth retardation. Both transgenic lines had no off-target mutagenesis in the investigated three genes ( Supplementary Fig. S7), though other mutations might occur in important genes by CRISPR/Cas9. On the other hand, we cannot deny the possibility of somaclonal mutation in GE#1-17 and GE#1-18. Few studies have been reported on the somaclonal mutation of sugi during cell culture and somatic embryogenesis, however, there are reports of somaclonal mutation or abnormal growth in other gymnosperms such as Picea [48][49][50] , Pinus 51 , and Larix 52,53 . Therefore, severe growth inhibition and lethality in GE#1-17 and GE#1-18 may be attributed to somaclonal mutation.
The sugi CjACOS5 gene was revealed to be essential for pollen production in this study. The Arabidopsis ACOS5 protein is an acyl-CoA synthetase that is estimated to be linked to the biosynthesis of sporopollenin in the exine of pollen 12 . The rice OsACOS12 protein is the ortholog of ACOS5 and catalyzes the condensation of oleic acid and CoA in vitro 14 . Although the physiological and biochemical functions of CjACOS5 were not clarified in this study, CjACOS5 was deduced to function as an acyl-CoA synthetase based on the structural analogy to ACOS5 and OsACOS12 (Fig. 1). Sporopollenin biosynthesis requires not only acyl-CoA synthetase, but also several other enzymes encoded by the MALE STERILITY 2 (MS2) 54 61,62 . The Arabidopsis ms2 mutant, the drl1/tkpr1 mutant, the lap3 mutant, and the double mutant of lap5 and lap6 (pksb and pksa) show complete male sterility. Similarly, the rice defective pollen wall (dpw, the ortholog of MS2) 63 , Ospks1 64 , Ospks2 65 , and Ostkpr1 66 mutants exhibit a male-sterile phenotype individually. The wheat TaNP1 gene (encoding a putative glucose-methanol-choline oxidoreductase) is an ortholog of OsNP1 that is involved in tapetum degeneration and pollen exine formation 67 . As the result of TaNP1 disruption by CRISPR/ Cas9, triple homozygous mutants of TaNP1 showed complete male sterility. The tomato SlSTR1 gene (a stamenspecific expressed strictosidine synthase) is an ortholog of Arabidopsis LAP3, and CRISPR/Cas9-mediated mutagenesis of SlSTR1 confers no viability on pollen 68 . Those male sterility-related gene orthologs of sugi probably can be utilized in addition to CjACOS5 to generate new varieties with male sterility by CRISPR/Cas9.
Our study established that the genome editing technology is available as a method for generating pollen-free sugi trees. So far, conventional crossbreeding programs using wild male-sterile sugi mutants have been useful for the generation of new pollen-free varieties of sugi. Because these crossed varieties can be freely handled, unlike the case of GMO trees, the new pollen-free sugi varieties that are the progeny of sugi plus trees and wild pollen-free mutants have been gradually planted in Japan 3 . In contrast, the CjACOS5-disrupted pollen-free sugi trees obtained in this study are GMO trees and, thus, are currently restricted for handling and planting by law. However, the genome editing technology, including CRISPR/Cas9, has several advantages compared with conventional tree breeding, as follows. First, if we could identify target genes involved in excellent traits, genome editing technology would be able to directly modify the target genes and generate their mutant trees in cases in which mutants of the target genes are not found in nature. Second, CRISPR/Cas9 can edit several target genes simultaneously and, thus, can modify multiple traits. Considering the combination of expected traits, such as no-pollen, fast growth, quality of wood, straightness, and stress tolerance, multiple values can be conferred to sugi breeding using multiplex gene editing technology. Last, the genome editing technology can not affect every gene other than the target genes, if off-target mutagenesis can be controlled fully. Therefore, the genome editing technology is difficult to perturb the genetic composition of well-established varieties unlike crossbreeding. In case of using genome-edited trees in forestry, transgene-free null segregants are actually needed. To generate null segregants, we have been crossing the CjACOS5-mutated sugi trees with non-transgenic sugi trees. By genetic recombination and segregation, in future, we expect to obtain null segregant sugi progenies without the transgenes such as NPTII and Cas9. Recently, several methods have been developed for the generation of nontransgenic genome-edited plants [69][70][71] . Consequently, it is expected that they will be utilized in combination with the advantage of genome editing technology and conventional breeding methods for tree breeding.

Methods
Vector construction. pZK_gYSA_FFCas9 vector 37 was digested with XbaI, the recognition site of which was located between the cauliflower mosaic virus 35S promoter and the NPTII gene, and blunted with Klenow fragment. About 2 kb of Zea mays ubiquitin promoter was inserted into the blunted pZK_gYSA_FFCas9 vector 72 . The resulting binary vector was designated as pBFGE1. Two synthetic DNAs (Supplementary Table S1) were annealed and ligated into the BbsI-digested pUC19_AtU6oligo vector 73 . A derivative vector was digested with I-SceI, and DNA fragments including the oligonucleotides were excised. The DNA fragments were inserted into I-SceI digested pBFGE1.
Plant materials and sugi transformation. The sugi (Cryptomeria japonica D. Don) embryonic cell line #13-8-2 was induced from an immature seed using previously described methods 74 . The #13-8-2 cells were transformed using A. tumefaciens GV3101 harboring the constructed pBFGE1 vector as described previously 75,76 . Selected kanamycin-resistant calli were regenerated through somatic embryogenesis 77 . Regenerated sugi plantlets were acclimatized in moist peat moss and Kanuma pumice, then transferred to moist vermiculite in a pot. They were cultivated under natural-day-length conditions in a closed glass greenhouse (phytotron) or a special netted-house. The temperature inside the phytotron was controlled by an air conditioner following the outdoor temperature. The special netted-house was a greenhouse with metal mesh windows that passed the outside air but impeded invasion by insects. www.nature.com/scientificreports/ DNA analysis. Orthologs of the CMFL003_A04 cDNA were searched using the blastn program in the latest NCBI nucleotide database 78,79 . The CjACOS5 cDNAs were cloned as described in Supplementary Information. A functional domain analysis of ACOS5-related protein sequences was performed on the InterPro website 80 . Genomic DNA was isolated from embryogenic cells and from leaves using the DNeasy Plant Mini Kit (Qiagen). The NPTII and CjACOS5 genes were confirmed using PCR (Supplementary Information). To elucidate the mutation induced by genome editing, a DNA region including the predicted DNA break site for the CRISPR/ Cas9 system was amplified using PCR (Supplementary Information). The amplified DNA fragments were cloned into the pBlueScript II vector and sequenced.
Induction and observation of male strobili. To induce flowering, 100 mg/L of a gibberellin A3 (GA 3 ) aqueous solution was sprayed on sugi plantlets. In the first experiment, GA 3 treatments were performed twice, on July 24 and August 3, 2017. Male strobili were observed in January, 2018. Pollen and the suspension of sliced male strobili were stained with Calberla's fuchsin staining solution, for microscopy 81 . In the second and third experiments, GA 3 treatments were performed in late July, 2018 and 2019, respectively. Thereafter, male strobili and pollen were analyzed in the following year in the same way as that described for first experiment.

Data availability
The data that support the findings of this study are available from the corresponding author upon reasonable request.